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ABSTRACT 

The  energetics  of  two  numerical  models  of  the  ocean 
circulation  are  studied  by  use  of  available  potential  energy. 
Available  potential  energy,  calculated  by  an  approximation 
equation  using  a  temperature  perturbation  field,  compares 
favorably  with  that  obtained  using  the  exact  expression 
using  a  height  perturbation  field. 

A  general  circulation  model  of  the  ocean  is  employed  to 
study  large-scale  energy  transformations,  and  a  mid-ocean 
dynamics  experiment  (mesoscale)  model  is  employed  to  study 
transformations  of  energy  by  stationary  features  and  by 
transient  motions.  Results  from  the  general  circulation 
model  are  consistent  with  those  of  Haney  (197^-).   Results 
from  the  mesoscale  model  indicate  that  the  transformation 
of  energy  is  from  available  potential  energy  to  kinetic 
energy  of  the  transient  eddy  field. 
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I.  INTRODUCTION 

A  "basic  part  of  any  study  of  an  atmospheric  or  oceanic 
circulation  system  is  a  study  of  the  behavior  of  the  energy 
in  the  system  and  its  transformation  from  one  form  of  energy 
to  another.   The  various  forms  of  energy  can  be  separated 
into  three  basic  forms  of  energy:  kinetic  energy,  potential 
energy  and  internal  energy.   Haurwitz  (1940  showed  that 
the  potential  energy  and  internal  energy,  in  a  vertical 
column,  are  proportional  under  hydrostatic  equilibrium  and 
thus  may  be  combined  into  a  single  parameter,  the  total 
potential  energy.   Conversions  of  energy  then  occur  between 
total  potential  energy  and  kinetic  energy. 

A.   AVAILABLE  POTENTIAL  ENERGY  IN  ATMOSPHERIC  CIRCULATION 

Lorenz  (19j?5)»  in  &  study  of  atmospheric  circulation, 

showed  that  since  total  potential  energy  was  three  to  four 

orders  of  magnitude  larger  than  kinetic  energy,  transformation 

of  energy  between  the  two  energy  forms  was  difficult  to  study. 

In  a  horizontally  stratified  atmosphere,  heating  of  the 

atmosphere  disturbs  the  stratification,  allowing  potential 

energy  to  be  converted  to  kinetic  energy.   However,  cooling 

of  the  atmosphere,  which  removes  total  potential  energy  from 

the  system,  also  disturbs  the  stratification  and  creates 

conditions  for  conversion  of  total  potential  energy  to  kinetic 
energy. 
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Lorenz  (1955)  defined  available  potential  energy  as  the 
difference  between  the  total  potential  energy  of  the  whole 
atmosphere  and  the  total  potential  energy  that  would  exist 
if  the  mass  of  the  atmosphere  were  redistributed  with 
conservation  of  potential  temperature  to  form  a  horizontally 
stable  atmosphere.   Then,  under  adiabatic  flow,  the  sum  of 
the  available  potential  energy  and  kinetic  energy  is 
conserved. 

B.   AVAILABLE  POTENTIAL  ENERGY  IN  OCEANIC  CIRCULATION 

The  use  of  the  concept  of  available  potential  energy 
can  be  carried  into  the  study  of  oceanic  circulations. 
Several  studies  of  the  energetics  of  the  ocean  circulation 
systems  have  been  based  on  available  potential  energy,  as 
in  Barrett  (1970,  and  Wright  (1972).  The  study  by  Wright 
dealt  with  changes  in  available  potential  energy  in  the 
ocean  general  circulation.  The  cooling  of  cold,  high  latitude 
water  and  the  surface  warming  of  warm  equatorial  water 
produces  available  potential  energy  by  increasing  the 
stratification  of  the  ocean  mass.   Conversely,  the  sinking 
of  the  cold  surface  water  in  the  north  tends  to  convert 
available  potential  energy  to  kinetic  energy  as  it  returns 
the  water  mass  to  a  more  stably  stratified  system.   Barrett 
(1971)  studied  the  available  potential  energy  in  the  transient 
motion  of  Gulf  Stream  rings.   Results  of  the  study  indicated 
that  the  decay  of  cold-core  rings  in  the  western  Sargasso 
Sea  region  reduced  the  available  potential  energy  of  the 
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region.   The  presence  of  decaying  cold-core  rings  in  the 
western  Sargasso  Sea  should  have  resulted  in  a  rising  of 
the  mean  thermocline  in  that  region.   Since  the  thermocline 
was  not  observed  to  rise  significantly  in  the  western 
Sargasso  Sea  region,  a  local  energy  source  must  have  been 
present  to  counteract  the  energy  loss. 


C.   AVAILABLE  POTENTIAL  ENERGY  IN  NUMERICAL  OCEAN 
CIRCULATION  MODELS 


This  study  applied  available  potential  energy  calculations 
to  two  experiments  with  a  numerical  model  of  the  ocean  general 
circulation;  one  simulating  large-scale,  quasi-steady  state 
motion  and  the  other  simulating  mesoscale  transient  motion. 
It  was  desired  to  calculate  available  potential  energy 
using  a  perturbation  height  field  compared  with  an  approximate 
calculation  using  a  temperature  perturbation  field.  An 
expression  for  the  balance  of  available  potential  energy  was 
obtained  from  this  approximation  equation.   The  importance 
of  the  terms  in  the  expression  for  the  balance  of  available 
potential  energy,  and  their  corresponding  physical  processes, 
could  then  be  examined  in  the  two  numerical  experiments. 
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II.  DESCRIPTION  OF  THE  NUMERICAL  MODELS 

The  first  experiment  utilized  the  general  circulation 
model  of  the  ocean,  as  described  in  part  A,  for  large-scale 
quasi-steady  state  motion.   The  Mid-Ocean  Dynamics  Experiment 
model  (Robinson,  et.al.,  1974,  MODE,  Hot  Line  News.  No.  57.), 
as  described  in  part  B,  was  utilized  for  the  transient  motion 
experiment. 

A.   GENERAL  CIRCULATION  MODEL  OF  THE  OCEAN 

The  model  represents  an  idealized  ocean  basin  with  a 
flat  bottom  and  regular  coast.   It  is  everywhere  4  km  deep, 
90°  longitude  wide  and  extends  from  51.25S  to  /+8.75N. 
The  grid  intervals  are  3°  of  longitude  and  2.5  of  latitude.  ■ 
There  are  six  vertical  levels,  placed  at  20,  100,  300,  700, 
1500,  3000  meters  below  the  undisturbed  surface.   The  model 
is  based  on  the  primitive  equations  which  are  integrated  in 
time  from  specified  initial  conditions.   The  equations  are 
based  on  the  hydrostatic  and  Boussinesq  approximations. 
The  ocean  is  thus  assumed  incompressible.   Salinity  is 
neglected  and  density  is  a  linear  function  of  temperature. 
The  equations  are 

dt   Pacosfd*       ac°sf  Oi3- 
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\\r  - 


dl         C0<*   £f  acosvp 


5e       -«  (3) 


£Z     ^      acos«f  L  2>a        <5f  J 


where  the  symbols  have  the  following  meanings: 

t  =  time 
*X  =  longitude 
f  =  latitude 
Z  =  height 

A  =  angular  speed  of  earth's  rotation 
f  =  coriolis  force  (  =  2  it  sin<^  ) 
a  =  radius  of  earth 
g  =  acceleration  of  gravity 
u  =  eastward  velocity  component 
v  =  northward  velocity  component 
w  =  vertical  velocity  component 
T  =  temperature 

T  =  constant  reference  temperature 
Q  =  density 


(4) 


(5) 


e=G(i-<0"-T.))  3  (6) 


p  =  density  of  water  at  reference  temperature 

eC  =  coefficient  of  thermal  expansion 

H  =  depth  of  ocean  basin 

e>lT)=  convective  adjustment  of  T  if  lapse 
rate  becomes  unstable. 

The  lateral  boundary  conditions  on  heat  and  momentum 

at  all  four  walls  of  the  basin  are  those  of  perfect  insulation 

and  zero  slip.   The  ocean  floor  is  thermally  insulated,  but 

free  slip  is  allowed.   Due  to  the  flat  bottom,  the  vertical 

velocity  is  zero  there.   The  circulation  is  driven  and 

controlled  by  the  wind  stress  and  heat  flux  at  the  ocean 

surface.   The  wind  stress  is  larger  in  the  southern  hemisphere 

and  meridional  wind  stress  is  neglected.   The  net  downward 

heat  flux  at  the  ocean  surface  is  expressed  as 

—•——■- 

where  T  is  the  ocean  surface  temperature,  T   is  an 

atmospheric  equilibrium  temperature  approximately  equal  to 

the  air  temperature  at  10  meters  and  Q~  is  a  coupling 

coefficient  determined  from  atmospheric  climatological 

data.   In  this  model,  T  is  the  temperature  at  the  top  level, 

20  meters  below  the  surface.   The  model  is  a  coarse-grid, 

high-viscosity  model  designed  to  show  large  scale  ocean 

circulation.   Complete  documentation  of  the  model  is  given 

in  Haney  (197*0. 
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B.   MID-OCEAN  DYNAMICS  EXPERIMENT  MODEL 

The  mid-ocean  dynamics  experiment,  MODE,  model  is  a 
modification  of  the  previously  described  model.   It  is 
designed  to  resolve  mesoscale  eddies  with  a  horizontal 
grid  resolution  of  ^-0  km  and  5  levels  in  the  vertical.   A 
2000  km  square  section  of  the  general  circulation  model  area 
is  represented  by  the  MODE  model,  located  along  the  western 
boundary  and  centered  at  35N.   Circulation  is  driven  and 
controlled  by  the  wind  stress  and  thermal  forcing,  the  same 
as  in  the  previous  model,  but  the  eddy  viscosity  and 
conductivity  were  reduced  to  achieve  mesoscale  eddies. 
Eastern  and  western  boundaries  are  thermally  insulated, 
no-slip  walls,  as  in  the  general  circulation  model.   Northern 
and  southern  boundary  conditions  are  taken  as  free  slip, 
zero  normal  velocity  and  a  vertical  temperature  distribution 
which  is  independent  of  longitude  and  time.   The  ocean 
floor  is  thermally  insulated  as  in  the  general  circulation 
model.   Initial  conditions  for  the  model  were  zero  velocity 
and  a  temperature  of  5  C  at  all  levels.   Twenty  years  of 

integration  was  performed  on  a  100  km  X  100  km  grid  using 

7   2 

a  horizontal  viscosity  of  A  =  2x10  cm  /sec  and  heat 

7   2  / 

diffusivity  of  A^  =  1x10  cm  /sec.   Ten  years  of  fine  grid 

6   2 

calculations  were  then  made  using  A  =  2x10  cm  /sec  and 

Ay  =  5x10   cm  /sec. 


16 


III.  DERIVATION  OF  AVAILABLE  POTENTIAL  ENERGY  EQUATIONS 

A.   DOMAIN  AVERAGED  AVAILABLE  POTENTIAL  ENERGY 

The  vertically  averaged  available  potential  energy 
contained  within  an  idealized  ocean  basin  may  be  calculated 
as  follows.   Assume  the  basin  to  be  flat  bottomed  at  a  depth 
H.  The  total  potential  energy  per  unit  volume  in  a  column 
in  the  basin  is  given  as 

p=  X  f  cgzdz  , 

where  Z  increases  upward  from  the  bottom  where  Z  =  -H.   All 
variations  in  H  are  neglected.   Neglecting  also  the  dependence 
of  density  on  pressure  and  salinity,  the  density  may  be 
written  as  a  linear  function  of  temperature,  equation  (6). 
The  total  potential  energy  per  unit  volume  in  the  column 
is  then  given  as  *s°  "H0)     -. 

By  defining  Z(T)  =  0  if  T  >  T(0)  and  Z(T)  =  -H  if  T  £  T(-H) , 
the  vertical  average  total  potential  energy  per  unit  volume 
in  the  column  may  be  written  as  _ 

zldr, 

where  Tmin  and  Tmax  represent  the  coldest  and  warmest 

temperatures  found  in  the  basin  respectively.   The  vertical 
average  total  potential  energy  averaged  over  the  horizontal 

domain,  denoted  by  (""*),  is  given  by 
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p  =  co^stawt   -  3oce0  t(-h)h*    ,     %*<?0  (   zz  ,T 

p 

where  the  bar  over  Z  denotes  an  average  over  an  isothermal 


surface.   The  average  potential  energy  is  a  minimum  when  Z 
is  a  minimum,  ie.  Z  is  equal  to  Z  everywhere.   Writing  the 
height  of  a  temperature  surface  as  a  mean  height,  Z,  and 
a  departure  from  the  mean,  Z1,  then  the  domain  average  total 
potential  energy  is  a  minimum  when  Z'  is  zero  everywhere.   In 
this  case,  Z  can  be  replaced  by  Z  .   The  minimum  domain 
averaged  total  potential  energy  can  then  be  written  as 

Tit  mi 

If  the  difference  between  the  average  bottom  temperature, 
T(-H) ,  and  the  average  bottom  temperature  in  the  state  of 
minimum  domain  averaged  total  potential  energy,  T(-H)min,  is 
assumed  small  and  thus  can  be  neglected,  then  the  domain 
average  available  potential  energy  can  be  written  as 


frj  ■&*'*)«> 


which  can  be  written  as 

J-r*dT.  (7) 

Thus,  the  available  potential  energy  is  proportional  to  the 
variance  of  height  of  an  isothermal  surface. 
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B.   APPROXIMATION  TO  THE  DOMAIN  AVERAGE  AVAILABLE  POTENTIAL 
ENERGY 

2 
The  variance  of  height  over  an  isothermal  surface,  Z1  , 

can  be  related  to  the  variance  of  temperature  over  a  horizontal 

surface,  assuming  the  ocean  is  stably  stratified.   When  the 

perturbation,  Z1,  from  the  mean  height  of  an  isothermal 

surface  is  greater  than  zero,  the  temperature  of  this  region 

is  less  than  the  mean  temperature  by  an  amount  T1.   Similarly, 

when  Z1  is  less  than  zero,  T1  is  greater  than  zero.   Therefore, 

changes  due  to  Z1  are  partially  compensated  by  changes  in  T1. 

The  height  of  a  temperature  surface  may  be  approximated  by 

Z(T)  %   1(f)  =Z(T-f)J 

which  shows  the  height  perturbation  can  be  approximated  by 

Z(T)  %    T'^I  , 

Y/hen  this  is  substituted  into  the  expression  for  the  domain 
average  available  potential  energy,  equation  (7),  the 
approximation  to  the  domain  average  available  potential 
energy  is  then 


) 


which  can  be  written  as 


(8) 
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C.   THE  BALANCE  OF  AVAILABLE  POTENTIAL  ENERGY 

Since  available  potential  energy  is  a  very  important 
source  for  kinetic  energy,  a  study  of  the  energetics  of  an 
ocean  circulation  should  include  a  study  of  production  and 
dissipation  of  available  potential  energy.   Under  adiabatic 
conditions  an  increase  in  kinetic  energy  is  compensated  by 
a  decrease  of  available  potential  energy.  When  both  energy 
forms  increase,  then  formation  of  available  potential  energy 
must  be  occurring  through  diabatic  processes.  A  study  of 
these  generation  and  dissipation  processes  and  their 
importance  may  be  made  through  use  of  the  equation  of  the 
balance  of  available  potential  energy.   The  equation  is 
derived  as  follows,  based  on  the  approximate  form  (8)  of 
available  potential  energy. 

From  equation  (8),  the  time  rate  of  change  of  available 
potential  energy  is 


hi   " 


-M(.ttJ  -H  \n) 

By  use  of  the  thermodynamic  equation  (5) >  equation  (9)  was 
rewritten  to  show  the  various  terms  responsible  for  the 
change  of  available  potential  energy.   The  term  Al'/«T'  J  can 


be  written  as 


Lfi.  T'1]  .  t'^I*  .  do) 

— •(  —  )  can  be  written  as 

MM-  UHY 
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The  local  rate  of  change  of  temperature  and  its  average  were 
obtained  from  the  thermodynamic  equation 

£>I  =  _V(VT)  _<^23  +  AHVZT  +  K  2l  \  &D.       (12) 
The  horizontal  average  of  (12)  gave 


^j  :^H     fc£i     s:ccr)  f  d3) 

£>Z      dZ       ^Z1        ' 

where  the  terms  -"^7*TVT)  and  JLj^T?  were  zero  since  advection 
and  diffusion  of  heat  vanish  at  the  boundaries  of  the  domain. 
Equation  (11)  was  then  written  using  (13)  as 


4(4f) ;  a  [_  S(wt)      n£j      4(t)1  m, 


which  can  be  rewritten  as 


where 


Equation  (12)  was  written  in  terms  of  an  average 
temperature  and  the  departure  from  the  average.   Equation  (13) 
was  then  subtracted  from  the  result  giving 

g=  -v(Vf)-v(Vf) -  yg -  ^  fjp^l^W ♦ ktf    ( . 6) 

Equation  (16)  was  then  substituted  into  the  right  hand  side 
of  equation  (10)  to  give 
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iffi  v\Hvn  ■  v}  f-^  -  ¥1  t  fgJM*} 


(17) 


since 


T  c*(v/t)  -  o. 


The  first  term  on  the  right  side  of  (17)  can  be  written  as 


tT-v(vt)-  SCwfjIs  -tVvt  -  t'tvv  -t'v/^x  -t'tsL*. 

|_  ^Z  J  3Z      2>Z 

The  term  -TV-VT  is  equal  to  zero  because  V  T  is  zero.  By 
use  of  the  continuity  equation  (k) ,    the  equation  reduces  to 


tT-  V'f^T)  -Mv»'t)1  -  -  t'w  df  .       (18) 
L         d"Z  J  "      3fc 

Also,  the  second  term  on  the  right  side  of  equation  (17) 
can  be  written  as 


TT-V-(VT')    -\CwVn      r    -  i-fwiS]    •  (19) 

The  last  term  of  equation  (17)  can  be  written  using  /the 
definition  of  Q  as 


Vf(AHvlT)'  +  (^/  +  5(T)'].-g 

Substitution  of  equations  (18),  (19),  (20)  into  equation  (17) 
and  using  the  resulting  equation  along  with  equation  (15) 
results  in  the  following  form  for  equation  (9),  the  balance 
of  available  potential  energy 
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(21) 


\n)  ISiJ   d 

The  balance  of  available  potential  energy  is  due  to 
five  terms  on  the  right  side  of  equation  (21).   The  first 
term  is  the  transformation  between  available  potential  energy 
and  kinetic  energy  due  to  the  correlation  between  temperature 
and  vertical  velocity.   The  second  term  is  a  triple 
correlation  term  due  to  the  convergence  of  the  vertical 
transport  of  temperature  perturbations.  The  third  term  is 
the  generation  of  available  potential  energy  due  to  a  spatial 
correlation  between  heating  and  perturbation  temperature. 
The  last  two  terms  represent  the  change  of  available  potential 
energy  due  to  a  change  in  the  mean  static  stability.   Lorenz 
(1955)  assumed  the  static  stability  was  constant  in  time  and 
also  neglected  the  triple  correlation  term  leaving  only  the 
first  and  third  terms  on  the  right  side  of  equation  (21). 
In  the  quasi  geostrophic  theory,  the  first  and  third  terms 
are  the  only  terms  to  appear  in  the  balance  of  available 
potential  energy.   For  the  general  circulation  model,  all 
five  terms  were  studied  and  for  the  mid-ocean  dynamics 
experiment  model  the  first  three  terms  were  studied. 
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IV.  FINITE  DIFFERENCE  EQUATIONS 

The  numerical  program  developed  for  calculations  of 
available  potential  energy  and  the  balance  of  available 
potential  energy  was  designed  for  the  general  circulation 
model  as  developed  by  Haney  (1974)«   Therefore,  the  space 
finite  differencing  scheme  and  the  boundary  conditions 
were  the  same  as  those  of  the  general  circulation  model, 
Haney  (1974).   A  staggered  grid  system  was  used  for  placement 
of  the  variables.   Figure  1  shows  the  variable  placement 
in  the  horizontal  and  figure  2  shows  the  variable  placement 
in  the  vertical.   Longitude  positions  are  designated  by 
the  subscript  i,  latitude  positions  by  the  subscript  j  and 
vertical  positions  by  k,  increasing  downward. 

Horizontal  averaging,  denoted  by  q  for  some  variable 
q. .  was  done  in  the  following  manner.   Since  the  lateral 
boundaries  were  located  at  integer  i,j  points,  the  mean 
value  of  a  term  to  be  averaged  was  calculated  at  the  i+-jf,  op- 
positions based  on  the  terms  value  at,  and  the  areas  of,  the 
grid  points  at  the  four  surrounding  i,j  positions.   These 

i 

mean  values  were  then  summed  with  equal  weights  to  determine 
the  domain  horizontal  average. 
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Figure  1.  Horizontal  Placement  of  Variables 
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Figure  2.  Vertical  Placement  of  Variables 
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A.   NUMERICAL  CALCULATION  OF  AVAILABLE  POTENTIAL  ENERGY  BY 
HEIGHT  PERTURBATION 

To  calculate  the  available  potential  energy  by  equation  (7) 

it  was  necessary  to  calculate  the  temperature  at  the  boundaries 

of  the  layers  represented  by  the  subscript  k.   T.  ..  represents 

the  vertical  mean  temperature  for  a  layer  of  thickness  DZ, 

v/here 


1  j.f2**  +4Z**J 


The  temperature  and  height  of  the  water  column  at  the 
boundaries  of  each  k  level  were  calculated  such  that  T..  .,, 
was  maintained  as  a  vertical  mean.   The  actual  temperatures 
9.  .,  i ,  were  calculated  using  the  following  equations 
6  ■  i  v  =    A-  Tea,  •  pZ,  -  0ija  '  P2»A 


where  KM  represents  the  total  number  of  levels  in  the  model. 
The  horizontal  average  of  the  vertical  mean  temperature,  T^  , 
was  calculated  as  explained  previously.  The  depth  of  T, 
was  determined  using  9.  .,  1  and  Z,  x   by  use  of  the   equation 

where  Z(T,  ).  ..  is  the  depth  at  which  the  temperature  9  is 
equal  to  the  horizontally  averaged  vertical  mean  temperature 
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in  the  layer  k.   The  horizontal  average  of  the  depth  of 
each  of  the  mean  temperature  surfaces  is  Z(T,),  ,  and 
represents  the  average  height  of  the  isotherms  between  9, 
and  9, + « .  The  departure  from  the  horizontal  average  was 
determined  as  the  difference  between  the  temperature  surface 
Z(T,) .  .,  and  its  horizontal  average  Z(T,  ),.   The  resulting 
finite  difference  analog  of  (7)  for  available  potential 
energy  was  of  the  form 


KK   -      _      


k--i 


B.   NUMERICAL  CALCULATION  OF  AVAILABLE  POTENTIAL  ENERGY  BY 
TEMPERATURE  PERTURBATION 

Equation  (8)  shows  the  analytical  forms  of  the  equation 
for  the  approximation  of  average  available  potential  energy 
using  the  temperature  perturbation,  T1,  instead  of  the  height 
perturbation,  Z1. 

The  horizontal  average  temperature,  T,  was  calculated  at 
each  k  level.   By  subtracting  this  value  from  the  temperature 
at  each  grid  position,  a  value  for  the  temperature  perturbation, 
T',  was  obtained  at  each  grid  position.   Using  these 
perturbation  temperatures,  the  term  iT"   was  calculated  at 
each  grid  position  and  averaged  to  the  half-integer  points 
(i+i,  j+i,  k+£). 

The  term  for  the  vertical  derivative  of  the  horizontal 
average  temperature  was  calculated  at  the  i,j  integer  points 
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as 

S1W»  ~    ^W,         z    -Zlw    » 

and  averaged  to  the  half-integer  points.  The  approximation 
to  the  mean  available  potential  energy  was  then  obtained  by 
combining  these  terms  and  summing  over  all  levels  as  indicated 

K,  

AW.fc&tMV***.  (22) 

t — .  / dT  \ 


C.   NUMERICAL  CALCULATION  OF  THE  BALANCE  OF  AVAILABLE 
POTENTIAL  ENERGY 


The  analytical  form  of  the  rate  of  change  of  available 

potential  energy  is  shown  by  equation  (9)  as 
_       —  ~*  -.       o 

6  a  _  3c<e 


-*1  Isi)  "H  l^/ 


»    "  «  ■  —  ■     ■  '  - ^    (9) 


The  factors  -^T1  and  /dU  in  (9),  were  calculated  as  in 
the  numerator  and  denominator  of  (22).   The  remaining  parts 
of  equation  (9)  were  calculated  by  use  of  the  right  side 
of  the  thermodynamic  equation  as  indicated  by  equations  (10) 

and  (11).   Equation  (10)  was  written  to  determine  ^r(^aT'  / 

St' 

by  calculating  T'  — -,     at  integer  values  of  i,j,k  and  then 

a  *■ 

averaging  to  obtain  its  value  at  half-integer  points. 
Equation  (11)  was  calculated  from 
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and  then  averaged  to  the  half-integer  grid  points.  The 
resulting  finite  difference  form  for  equation  (9)>  the 
balance  of  domain  average  available  potential  energy,  is 
then 

S5-.  fa&f  f-L  h[uf)  -  (kr.V.(il)  >L-  (33) 
at       h  L-  LfamtA*    /      /St\j    azlaw  J   *'4 


D.   NUMERICAL  CALCULATIONS  OF  THE  TERMS  FOR  THE  BALANCE 
OF  AVAILABLE  POTENTIAL  ENERGY 


As  indicated  by  equation  (21),  the  balance  of  mean 
available  potential  energy  contains  five  terms  which  are 
obtained  by  use  of  the  thermodynamic  equation  (5)» 
Development  of  equation  (21)  showed  that  the  five  terms 
could  be  written  as  below 


»  L  &  J 


(18) 


a/.L©)-.   IZ^^i  (24) 
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The  terras  in  the  brackets  on  the  right  side  of  the  above 
five  equations  were  calculated  using  the  identical  finite 
difference  method  as  used  in  the  general  circulation  model 
to  calculate  the  terms  of  the  thermodynamic  equation.   The 
details  are  in  Haney  (1974). 


31 


V.   RESULTS 

To  make  an  accurate  comparison  of  the  two  expressions 
for  domain  average  available  potential  energy,  equations  (7) 
and  (8),  and  to  determine  the  importance  of  the  terms 
occurring  in  the  balance  of  available  potential  energy, 
equation  (21),  the  general  circulation  model,  as  described 
in  chapter  II,  was  run  for  large-scale  motion.  A  description 
of  the  conditions  and  the  results  obtained  are  given  below. 

A.   GENERAL  CIRCULATION  MODEL  OF  THE  OCEAN 

To  establish  quasi-steady  state  conditions,  the  model, 
Haney  (1974),  was  run  in  two  phases.  For  the  first  phase, 
the  linearized  form  of  the  vertical  shear  current  was  used. 
The  vertical  shear  current  equations  were  obtained  from  the 
primitive  equation  form  of  the  horizontal  equations  of  motion, 
(1)  and  (2),  by  subtracting  the  vertical  mean  of  equation  (1) 
from  equation  (1)  and  by  subtracting  the  vertical  mean  of 
equation  (2)  from  equation  (2).   The  vertical  mean  current  is 
assumed  constant  in  time.   The  temperature  is  determined  from 
the  complete  thermal  equation  ( 5) •   The  first  phase  integration 
was  for  approximately  100  years.   For  the  second  phase,  the 
complete  nonlinear  equations  were  used.   The  second  phase 
integration  was  for  25  years.   A  more  complete  description  of 
phase  1  and  phase  2  may  be  found  in  Haney  (1974).   Twenty- four 
sets  of  model  history  data  was  taken  at  this  point  at  intervals 
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of  25  days,  representing  a  total  of  600  days.   The  results 
showed  the  motion  to  be  not  strictly  steady-state.   As  shown 
by  Haney  (197A-),  slight  warming  was  occurring  near  the  bottom 
after  the  integration  period  indicating  a  continuing  downward 
flux  of  heat.   Figure  3  shows  the  temperature  at  20  meters 
below  the  surface,  averaged  over  the  last  year  of  integration 
for  phase  2.   Figure  4  shows  the  vertical  velocity  at  60  meters 
below  the  surface  averaged  over  the  last  year  of  the  second 
phase.   Comparing  figures  3  s.nd  k   shows  two  large  gyres 
centered  on  each  side  of  the  equator.   In  these  gyres,  warm 
water  is  being  forced  down.   Near  the  equator  and  along  the 
northern  and  southern  boundaries,  cold  water  is  being  forced 
upward.   Thus,  available  potential  energy  is  being  generated 
by  these  motions.  j,  
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Figure  3.  Temperature  (°C)  at  20  m,  averaged  over 
tne'last  year  of  the  second  phase. 
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Figure  ^.  Vertical  velocity  (cm  day  )  at  60  m, 
averaged  over  the  last  year  of  the 
second  phase.  Upward  motion  is 
shaded. 


1 .   Results  of  Available  Potential  Energy  Calculations 

Comparison  of  calculations  of  available  potential 
energy  as  determined  by  equation  (8)  and  the  approximate 
form  (9)  showed  that  the  available  potential  energy  as 
calculated  by  the  height  perturbation  to  be  slightly  larger, 
(1383  ergs/cm-5) ,  than  the  available  potential  energy  as 

•z 

calculated  by  the  temperature  perturbation,  (1232  ergs/cnr ) , 
figure  5.   This  discrepancy  is  perhaps  due  to  Tmax  and  Tmin 
being  not  uniquely  defined  in  Z  coordinates.   Fluctuations 
in  available  potential  energy -due  to  a  minor  change  in  the 
model  near  the  end  of  the  600  day  period  of  record  were 
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Figure  5.  APE  by  iZ,2(»)  and  by  iT,2(.)  for  the 
General  Circulation  Model. 
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simulated  by  both  forms  of  available  potential  .energy 
calculations.   Figure  5  indicates  equation  (8)  calculated 
available  potential  energy  to  a  sufficient  degree  of  accuracy 
to  be  considered  a  good  approximation  equation  for  available 
potential  energy  in  this  model. 

2.   Results  of  Calculations  of  the  Balance  of 
Available  Potential  Energy 

From  equation  (21),  the  balance  of  available  potential 

energy  is  represented  by  the  sum  of  two  integrals, 

&A  -   part  1  -  part  2 
St 

where  0  


parti  --  3f<e»/T-T,W-J_.iL  (*i:*\*  -i-LSl 

»  JA        [S|»u   I    pf\e0cj 

"    i.    wl    2.        W3. 


(26) 


par+a  -    ^S 


According  to  their  definitions,  Q1  and  Q  contain  the  convective 
adjustment  of  temperature,  d„(T)',  and  % „(T)  respectively. 
In  the  calculations  below,  these  convective  adjustment  terms 
were  omitted  due  to  the  difficulty  of  recovering  this 
irreversible  process  from  the  history  records. 

Part  1  represents  changes  in  available  potential 
energy  resulting  from  changes  in  the  variance  of  the  temperature 
perturbation  on  a  horizontal  surface.   Heating  of  warmer 
regions  and  cooling  of  colder  regions  at  the  same  depth 
increases  the  horizontal  temperature  variance  and  increases 
available  potential  energy.   Part  2  represents  changes  in 
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available  potential  energy  due  to  changes  in  mean  static 

stability.   Heating  of  upper  levels  and  cooling  of  lower 

levels  increases  the  static  stability  and  decreases  available 

potential  energy.   The  magnitude  of  these  integrals  are 

shov/n  by  figure  6.   Variations  in  these  two  parts  after 

day  ZfOO  of  the  history  record  were  caused  by  the  fact  that 

the  convective  adjustment,  $    (T) ,  was  changed  to  double 

precision  accuracy  in  the  model. 

The  processes  represented  by  part  1  are  serving  to 

decrease  the  total  available  potential  energy  in  this  model, 

indicating  the  model,  apart  from  &  (T) ■ ,  is  working  to 

c 

reduce  temperature  perturbations.   The  magnitude  of  the 
three  terms  comprising  part  1,  equation  (26),  are  shown 
in  figure  7,  where  all  three  terms  are  of  the  same  order  of 
magnitude.   The  first  term  of  part  1  shows  the  transformation 
of  energy  between  available  potential  energy  and  kinetic 
energy  due  to  a  correlation  between  the  temperature 
perturbation  and  vertical  velocity.   Results  indicate  that 
in  the  domain  average,  energy  is  being  transformed  from 
kinetic  energy  to  available  potential  energy.   This  is 
consistent  with  the  results  of  Haney  (197^)>  as  discussed 
above,  where  a  transformation  of  kinetic  energy  of  vertical 
shear  motion  to  potential  energy  was  occurring.   Term  2 
represents  a  triple  correlation  due  to  convergence  of  the 
vertical  flux  of  temperature  fluctuations.   Being  a  triple 
correlation  of  perturbation  quantities,  Lorenz  (1955) 
assumed  this  term  to  be  negligible  in  the  atmosphere.   In 
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Figure   6.   The  Values  of   (Part1-Part2,   •)    and 
the  Values  of  Parti (■)    and  Part2(*) 
of  Eq.    (26)    vs.    Time. 
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Figure  7.  Contribution  by  Term1(»), 
Term2("),  Term3(#)  to  the 
Balance  of  APE,  Eq.  (26). 
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this  general  circulation  model,  term  2  is  smaller  than  all 
of  the  other  terms  by  only  a  factor  of  2  to  3  and  thus  it  cannot 
be  neglected.   This  term  arises  because  the  static  stability, 
^T/^Zj  is  not  constant  in  Z.   If  it  were,  the  term  would 
vanish.   The  term  would  also  vanish  if  the  triple  correlation 
were  small.  One  reason  it  was  not  small  is  because  of  the 
forced  nature  of  the  vertical  velocity,  w,  and  temperature 
fields  near  the  boundaries.  With  heat  flux  through  the 
boundaries,  neither  the  vertical  velocity  nor  the  temperature 
gradient  in  the  horizontal  would  have  been  as  large  as 
occurred.  The  third  term  of  part  1  shows  generation  of 
available  potential  energy  due  to  a  spatial  correlation 
between  heating  and  perturbation  temperature.  This  term 

7/as  calculated  without  the  contribution  from  S_(T)', 

c 

present  in  equation  (20).  The  negative  result  of  term  3 
indicates  a  loss  of  available  potential  energy  due  to  the 
non-convective  heating.  The  magnitude  of  this  term  indicates 
the  importance  of  the  thermal  diffusion  terms  in  this 
oceanic  general  circulation  model. 

The  processes  represented  by  the  terms  of  part  2, 
figure  8,  were  of  the  same  order  of  magnitude,  and  tend  to 
increase  the  available  potential  energy  in  the  model.  The 
first  term  of  part  2,  term  l±t   is  due  to  the  vertical  flux 
of  mean  temperature  and  the  second  term,  term  5>  is  due  to 

mean  heating.  As  in  term  3>  the  cdnvective  adjustment 

term,  3  (t) ,  was  absent  from  the  mean  heating  term.  $  (T) 
c  c 

tends  to  increase  the  mean  static  stability,  £>T/^Z,  and  to 
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Figure  8.  Contribution  by  TernUf(») 
and  Term5(*)  to  the 
Balance  of  APE,  Eq.  (26) 


M 


decrease  available  potential  energy.   It's  inclusion  would 
allow  terms  Zf  and  5  to  sum  to  a  negligible  value.   Lorenz 
(1955)  assumed  part  2  was  zero  for  steady-state  atmospheric 
motion.   Under  conditions  of  seasonal  radiation  heating, 
the  vertical  flux  of  mean  temperature,  term  Zf,  may  not  be 
in  phase  with  the  mean  heating  that  occurs,  term  5>  and 
thus  would  not  cancel.   The  general  circulation- model  was 
not  designed  to  show  seasonal  radiation  heating. 

The  energy  flow  diagram  depicting  terms  1  -  5  in 
equation  (26)  is  shown  in  figure  9»  Terms  1,4  and  5  were 
serving  to  increase  available  potential  energy  by  an  amount 
of  1.705x10  J   ergs/cm  sec.  Terms  2  and  3  were  serving  to 
decrease  available  potential  energy  by  an  amount  of  .908x10  J 
ergs/cur  sec.   The  available  potential  energy  during  the 
record  period  was  observed  to  be  decreasing,  figure  5,  with 
an  average  value  of  1232.  ergs/cm  .   The  discrepancy  between 
the  terms  serving  to  increase  available  potential  energy, 
part  1  -  part  2)0,  and  the  observed  decrease  in  available 
potential  energy,  figure  5,  was  determined  to  be  due  to  not 
including  the  o.(T)  term  in  the  balance  of  available  potential 
energy  calculations.   The  value  of  the  term  required  to 
balance  the  energy  flow  diagram,  as  in  steady-state  conditions, 
is  indicated  by  term  6  and  shown  by  the  dashed  line  in 
figure  9  to  be  approximately  .797x10  J   ergs/cnr  sec.  These 
results  indicated  convective  adjustment  was  an  important 
process  in  the  model,  as  was  expected  since  the  model 
contained  a  large  region  of  near  isothermal  water  at  high 
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Figure  9»  Energy  Flov/  Diagram  for  the 
General  Circulation  Model. 
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latitudes  where  convective  adjustment  would  frequently 
occur.  . The  convective  overturning  plays  no  role  in  the 
kinetic  energy  budget  and  therefore  the  arrow  for  term  6 
shows  this  effect  to  be  an  external  loss  to  the  system. 

B.   MID-OCEAN  DYNAMICS  EXPERIMENT  MODEL 

Figure  10  below  shows  the  instantaneous  streamf unction 
for  day  10080  during  the  seventh  year  of  fine  grid 
calculations.  Twenty- four  sets  of  data  for  the  experiment 
were  taken  from  day  10050  to  10110  at  2?   day  intervals. 


Figure  10.  Instantaneous  streamfunction  for 
day  10080  of  the  Mid-Ocean 
Dynamics  Experiment  Model. 
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To  determine  the  available  potential  energy  as  a  result 
of  transient  motion  and  standing  motion,  equation  (8)  was 
written  with  temperature  in  terms  of  a  time  average 

temperature,  £t] ,  for  standing  motion,  and  a  departure  from 

* 

the  time  average,  T  ,  for  transient  motion.   Since  available 

potential  energy  is  determined  from  the  temperature 
perturbation,  the  horizontal  average  of  the  time  averaged 

temperature  was  determined.   Then  the  perturbation  temperature 

*  i 

of  the  standing  motion,  [T] •  and  transient  motion,  T   ,  was 

determined  as  below 

where  \JI^   is  the  horizontal  averaged  standing  temperature 

and  T  is  the  horizontal  averaged  transient  temperature. 

r  *t 

The  term  \.T  j  '  is  the  standing  eddy  temperature  and  T   is 

the  transient  eddy  temperature.   The  eddy  temperature  is 
then  the  sum  of  the  last  two  terms  of  equation  (2.7). 
Substitution  of  the  eddy  temperatures  into  equation  (8) 
resulted  in  the  following  form  for  domain  average  available 
potential  energy, 


h5 


The  time  averaged  available  potential  energy  was  written 
as 


where  the  term  JjKTl'T  ]  vanishes.   In  this  model,  the 
static  stability  factor  Lei  T/c>Z]  was  assumed  constant  in 
time  and  written  as  "£  .   Equation  (29)  shows  time  average 
available  potential  energy  as  the  sum  of  a  part  due  to  the 
standing  field  and  a  part  due  to  the  transient  motion. 

For  the  study  of  the  balance  of  available  potential 
energy  in  the  mid-ocean  dynamics  experiment  model,  only 
the  first  three  terms  of  equation  (21)  were  considered 
important  for  further  study.   Because  lateral  boundary 
conditions  allowed  heat  to  flow  in  or  out  of  the  model,  it 
was  reasoned  that  convective  adjustment  would  not  play  an 
important  role  in  the  model.   With  little  convective 
adjustment  and  steady  forcing,  it  was  therefore  assumed 
that  the  long-time  average  of  the  sum  of  terms  ^  and  5» 
part  2  of  equation  (26),  would  be  negligible.   Calculations 
of  the  balance  of  total  available  potential  energy,  below, 
showed  the  terms  of  part  2  summed  to  approximately  zero 
for  each  set  of  data  as  well  as  in  a  time  average  so  the 
assumption  was  justified.   The  study  therefore  concentrated 
on  transformations  of  energy  due  to  iT1  ,   From  the  first 
three  terms  of  equation  (21),  the  balance  of  available 
potential  energy  was  approximated  as 


^6 


I  .1 

(30) 


1.         2.  3. 

Equation  (30)  was  written  in  terms  of  a  part  due  to  transient 

*  r-  -• 

motion,  q  ,  and  a  part  due  to  standing  motion,  ^qj ,  for  a 

quantity  q  and  then  time  averaged  as 


UrJ  "  "TT7  1  **A2  2 


1a.       1b.  2a.        2b. 


2c.        2d.         3a.       3b. 
where  the  quantities  containing  only  one  term  due  to  transient 
motion  were  zero  in  the  long-time  average.   The  values  of 
terms  2c  and  2d  were  not  calculated  in  this  study  but  their 
sum  could  be  inferred  from  terms  2,  2a  and  2b.   Since  heat 
diffusivity  was  allowed  across  the  boundaries,  the  definition 
of  Q/Cftc  in  (25)  was  changed  to  the  following 
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o 

since  A^V  T  was  not  equal  to  zero  in  the  MODE  model. 

1 .   Results  of  Available  Potential  Energy  Calculations 


Figure  1 1  shows  the  total  available  potential  energy 
as  determined  from  equation  (28)  to  be  slowly  increasing 
during  the  history  record.   The  first  term  on  the  right  side 
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Figure  11.  APE  vs.  Time  for  the  MODE  Model 
(Total  Motion) 
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of  equation  (28),  due  to  standing  motion  was  constant  in  time 
and  had  a  value  equal  to  116.65  ergs/cm  ,  shov/n  in  figure  11. 
This  was  an  order  of  magnitude  smaller  than  the  available 
potential  energy  of  the  general  circulation  model  and  it 
points  out  the  importance  of  the  horizontal  diffusion  of 
heat  across  the  northern  and  southern  boundaries  in  preventing 
large  horizontal  gradients  of  temperature.   The.  third  term 
of  equation  (28) ,  due  to  transient  motion,  is  shown  in 
figure  12.   This  term  appears  to  have  a  period  of  approximately 
100  days,  but  is  too  small  in  amplitude  to  significantly 
affect  the  total  available  potential  energy.   Since  the  left 
side  of  equation  (28)  was  known  from  figure  11,  and  two 
terms  on  the  right  side  of  equation  (28)  were  known,  the 
middle  term  on  the  right  side  of  equation  (28)  was  found  to 
be  primarily  responsible  for  the  slow  increase  of  available 
potential  energy  seen  in  figure  1 1 .   The  time  average  of 
the  total  available  potential  energy,  corresponding  to  the 
left  side  of  equation  (29)  was  117.3  ergs/cnr.   The  first 
term  on  the  right  side  of  equation  (29),  the  available 
potential  energy  due  to  standing  motion,  is  the  same  as  the 
first  term  on  the  right  side  of  equation  (2.8),    equal,  to 
116.65  ergs/cm*.   Thus  the  time  average  available  potential 
energy  due  to  transient  motion,  the  last  term  of  equation  (29), 
had  a  value  of  .66  ergs/cnr  and  was  only  a  small  part  of  the 
time  average  total  available  potential  energy.   Being  so 
small,  the  available  potential  energy  due  to  transient 
motion  would  probably  not  be  noticed  on  a  temperature  map 


k9 


ergs 


cnr  sec 


0.0.  _ 


1.0 


30.0 

DAY  of  RECORD 


60.0 


Figure  12.  APE  by  £T»2  vs.  Time  for  the  MODS  Model 
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from  the  model.   If  the  variations  in  kinetic  energy  were  of 
this  order  of  magnitude,  it  could  be  important  in  the  kinetic 
energy  balance  which  would  be  subject  to  future  investigation. 

2.   Results  of  the  Balance  of  Available  Potential 
Energy  Calculations 

By  the  use  of  flow  diagrams,  the  magnitude  of  the 
terms  responsible  for  the  balance  of  available  potential 
energy  can  be  easily  seen.   Figure  13  shows  the  balance  of 
the  time  average  total  available  potential  energy, 
corresponding  to  equation  (30) •   This  figure  indicates  energy 
was  being  transformed  from  kinetic  energy  to  available 
potential  energy  by  term  1  of  equation  (30)  at  a  rate  faster 
than  it  was  being  removed  by  terms  2  and  3»  The  average 

rate  of  increase  of  total  available  potential  energy,  as 

—5       3 

shown  in  the  diagram,  was  approximately  .0^-6x10   ergs/cm  sec. 

From  the  observed  change  in  available  potential  energy, 

figure  1 1 ,  there  was  an  average  rate  of  available  potential 

—5        3 
energy  change  during  the  60  days  of  .03x10  J   ergs/cm  sec. 

Figure  1^-  shows  the  value  of  the  energy  transformations 

due  to  standing  motion,  corresponding  to  terms  1a,  2a,  3a  of 

equation  (31 )>  and  figure  15  shows  the  value  of  the  energy 

transformations  due  to  transient  motion,   corresponding  to 

terms  1b,  2b,  3b  of  equation  (31  )•   The  average  rate  of 

increase  of  available  potential  energy  by  the  standing  field, 

-5 
as  shown  in  the  diagram,  was  approximately  .075x10  ' 

ergs/cm  sec.   The  average  rate  of  decrease  by  the  transient 

field,  as  shown  in  the  diagram,  was  .035x10   ergs/cnr  sec. 
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Figure  13.  Energy  Flow  Diagram  for  the  MODE  Model 
(Total  Motion  Eq.  30) 
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Figure  1*f.  Energy  Flow  Diagram  for  the  MODE  Model 
(Standing  Motion,  Eq.  31) 
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Figure  13.  Energy  Flow  Diagram,  for  the  MODE  Model 
(Transient  Motion,  Eq.  30 
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Instead  of  such  changes  in  the  available  potential  energy 
fields,  this  perhaps  simply  indicates  that  a  conversion  of 
available  potential  energy  from  the  standing  to  the  transient 
fields  was  occurring.  The  conversion  term  for  available 
potential  energy  between  the  standing  and  transient  fields 
was  not  calculated  because  it  also  contained  terms  of 
standing  and  transient  motion  as  does  terms  2c  and  2d.   The 

conversion  term  must  have  had  a  value  between  .035  to  .075 

-5/3 
x10  v  ergs/cnr  sec.   Summing  of  the  net  increases  by  the 

standing  field  and  the  transient  field  gave  a  net  increase 

—  5        -z 
for  the  total  of  .040x10  y  ergs/cnr  sec,  which  compares 

with  the  value  of  .0^6x1  O*"'7  ergs/cnr  sec  seen  in  figure  13. 

The  sum  of  terms  3a,  figure  1Zf,  and  3"b>  figure  15,  was  equal 

to  term  3  of  figure  13  as  it  should  be,  since  in  the  time 

average  term  3  is  composed  of  only  term  3a  and  3t>.  Similarly, 

term  2a,  figure  14,  and  2b,  figure  15,  sum  to  equal  term  2 

of  figure  13.  However,  from  equation  (31)  term  2  is  composed 

of  four  terms  in  the  time  average.  Therefore,  terms  2c  and 

2d  must  be  negligible  or  sum  to  a  negligible  value. 

The  sum  of  la  and  1b  is  less  than  term  1  by  an 

—5        3 

amount  of  .006x10  J   ergs/cnr  sec.   Some  round-off  errors 

may  have  occurred  to  cause  the  small  discrepancy.   Comparing 
figure  13  and  figure  14,  the  total  energy  being  transformed 
from  kinetic  energy  to  available  potential  energy,  term  1 , 

was  almost  totally  by  the  standing  field,  term  1a.  Since 
term  1b,  figure  15,  was  small  and  positive,  the  source  of 
kinetic  energy  for  the  transient  motion  was  not  available 
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potential  energy  during  the  60  day  period  of  the  MODS  model. 
Figure  16  shows  the  total  transformation  from  kinetic  energy 
to  available  potential  energy,  term  1 ,  was  undergoing  large 
fluctuations.   At  some  times  the  transient  eddies  did  draw 
energy  from  available  potential  energy,  since  1b,  figure  17, 
was  at  times  negative,  but  the  energy  flow  in  either  direction 
was  small.   These  results  may  indicate  that  the- transient 
eddies  were  not  baroclinic  for  the  60  day  period.   The 
amplitude  of  the  eddies  of  the  MODE  model  were  known  to  be 
large.   With  large  amplitude  eddies,  the  dissipative  forces 
were  expected  to  be  large,  indicating  the  transformation  of 
energy  from  the  source  of  energy  for  the  kinetic  energy  of 
the  eddies  was  large.  Therefore  the  primary  source  of  energy 
for  the  transient  eddies  was  not  available  potential  energy 
during  the  60  day  history  record. 

It  was  suspected  that  errors  may  have  occurred  due 
to  the  2-*f  day  rate  of  sampling  the  data.   Figure  16  shows  the 
three  terms  of  equation  (30)  as  a  function  of  time.   Term  1 
and  2  fluctuate  to  a  greater  extent  than  term  3»   Since 
both  term  1  and  2  are  correlated  with  vertical  velocity,  w, 
the  possibility  of  aliasing  is  further  supported.   Fluctuations 
in  term  lb,  figure  17,  and  2b,  figure  18,  also  lends  support 
to  this  possibility  due  to  their  fluctuations.   Neither 
term  3  of  figure  16  nor  term  $b   of  figure  19  showed  a  large 

degree  of  fluctuation  and  neither  of  these  terms  were 
correlated  with  vertical  velocity.   Therefore,  the  above 
energetics  were  recalculated  using  a  120  day  history  record 
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Figure   16.    Contribution  by  Terni1x1Cr  (•) , 
Term2x106(«),    Term3x105(#) ,    to   the 
Balance   of  APE,    Eq.    (30). 
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Figure  17.  Contribution  by  Termlb  x10' 
to  the  Balance  of  APE. 
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which  began  on  day  10050  and  ended  on  day  101 70.   No  essential 
changes  occurred  in  any  of  the  terms  except  term  1b,  which 
was  still  small,  but  negative  for  the  120  day  average.   Since 
a  120  day  average  would  probable  be  more  accurate  than  a  60 
day  average,  the  eddies  were  felt  to  be  baroclinic. 
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VI.  CONCLUSION 

The  results  obtained  from  this  study  indicated  that 
both  the  general  circulation  model  and  the  MODE  model  were 
mechanically  driven  by  the  wind  stress  and  not  thermally 
driven  to  a  great  extent.   The  below  diagram  shows  the 
transformation  of  energy  from  its  generation  by  the  wind 
stress  to  its  dissipation  through  transient  eddies  as 
indicated  by  the  results  from  the  MODE  model. 
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The  atmospheric  wind  circulation  acts  on  the  ocean  surface 
creating  large  scale  standing  motion.   Some  of  the  kinetic 
energy  of  the  standing  motion  is  transformed  to  standing 
available  potential  energy  of  the  ocean.  Approximately 
half  of  the  input  of  energy  to  standing  available  potential 
energy  is  lost  through  non-convective  heating  by  standing 
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motion.   About  one-fourth  of  standing  available  potential 
energy  is  transformed  to  transient  available  potential  energy, 
where  a  large  portion  is  lost  through  transient  non-convective 
heating.   A  small  portion,  in  the  transient  available 
potential  energy  balance,  is  transformed  to  kinetic  energy 
of  transient  motion  where  it  is  lost  through  dissipation. 
V/hether  the  transient  available  potential  energy  is  the 
primary  source  for  the  transient  kinetic  energy  must  be 
determined  by  a  study  of  the  transient  kinetic  energy 
balance. 

In  the  atmosphere,  forcing  is  by  the  heating  and  the 
energy  transformations  are  from  available  potential  energy 
of  standing  motion  to  transient  available  potential  energy, 
then  to  transient  kinetic  energy. 

The  main  conclusion  is  that  baroclinic  instability  is 
a  possibility  for  maintaining  the  MODE  eddies.   The  general 
circulation  model  showed  energy  transformations  from  eddy 
kinetic  energy  to  eddy  available  potential  energy  for  large- 
scale  quasi-steady  state  motion.   The  MODE  model  showed  the 
energy  transformations  to  be  opposite,  from  eddy  available 
potential  energy  to  eddy  kinetic  energy,  for  transient 
motion.   Therefore,  during  periods  of  baroclinic  instability, 
eddy  available  potential  energy  is  aiding  in  the  maintenance 
of  the  eddies. 
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